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We find the perturbation spectrum of a family of spherically symmetric and continuously self- 
similar (CSS) exact solutions that appear to be relevant for the critical collapse of scalar field matter 
in 2+1 spacetime dimensions. The rate of exponential growth of the unstable perturbation yields the 
critical exponent. Our results are compared to the numerical simulations of Pretorius and Choptuik 
and are inconclusive: We find a CSS solution with exactly one unstable mode, which suggests that 
it may be the critical solution, but another CSS solution which has three unstable modes fits the 
numerically found critical solution better. 

PACS numbers: 04.20.-q, 04.25.-g, 04.40.-b 



I. INTRODUCTION 



Critical gravitational collapse, as first found by Choptuik JjJ, has been studied in many systems ||. For the most 
part, these studies have been numerical: generally the critical solution has not been found in closed form. However, 
gravitation in 2+1 dimensions is a much simpler system than in 3+1 dimensions. Thus one might hope to find 
a closed form solution of 2+1 critical collapse. Pretorius and Choptuik H performed numerical simulations of 2+1 
• critical gravitational collapse with a massless, minimally coupled scalar field and a cosmological constant. A numerical 
treatment of this system has also been done by Husain and Olivier. j| One of us || found a closed form continuously 
. self-similar (CSS) solution of the 2+1 Einstein-scalar equations that agrees with the work of Ref. ||. Other closed 
form solutions of these equations have also been treated^, [jj and the approach to the singularity in this system has 
O [ been analyzed. || 

A critical solution, when perturbed, should have exactly one unstable mode, which grows as e kT for some constant 
^ [ k, where T is a coordinate such that d/dT is the homothetic vector of the CSS critical solution. In the collapse of a 
one parameter family of initial data, a quantity Q with dimension (length)* obeys a scaling relation Q cx \p — p*\ ls 
where p is the parameter and p* is its critical value. The quantities k and 7 are related by 7 = 1/fc. Thus by treating 
• . perturbations of a critical solution, one could find k and compare to the numerical value of 7 found in simulations of 
near critical collapse. Such a perturbation treatment was begun in || but was not completed due to questions about 
appropriate boundary conditions to impose on the perturbations. 

In this paper, we complete the perturbation treatment begun in §. Throughout we use double null coordinates 
rather than the Bondi coordinates used in || . This tends to clarify the issue of boundary conditions for the perturba- 
tions. Section II presents the field equations and background solution in double null coordinates. The perturbations 
are treated in section III. Conclusions are presented in section IV. 

II. FIELD EQUATIONS AND BACKGROUND SOLUTION 

The simulations of Ref. [g[ used the Einstein-scalar equation with cosmological constant. However, in the approach 
to the critical solution, the cosmological term becomes negligible. Therefore we approximate the Einstein equations 
as 

R ab = 47rV a 0V 6 0. (1) 
We consider an axisymmetric 2+1 metric in double null coordinates, which takes the form 

ds 2 = 2e 2A dudw + B 2 d9 2 , (2) 
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where A and B are functions of u and w. We choose u at the origin to be proper time from the singularity, but place 
no restrictions on w. For a metric of the form (H), the Einstcin-scalar equations (Qj) become 
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(plus two additional components of the field equations which are redundant). The quantity A is fixed at the origin 
by the requirement that u be proper time there. 

To present the background solution, it is helpful to introduce coordinates T and y given by — u = e~ T and 
y = w n I 'y 1 —u where n is a constant. There is a family of regular CSS solutions parameterized by a positive integer n 
(the q of Ref. @), with 
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Here, the constant c is given by c — ±^J(2n — l)/87m and a bar denotes a background quantity. The solution with 
n — 4 appears to be a good fit to the critical solution found by Pretorius and Choptuik ^ in numerical collapse 
simulations. 

The spacetime is CSS with homothetic vector d/dT. The log-scale coordinate T used here is the same as that 
defined in Ref. B. The self-similarity coordinate y defined here is related to the coordinate R defined in Ref. (5) by 
y 2 = 1 — 2R. However, this simple relation between y and R does not hold when the spacetime is perturbed. Note 
that in these coordinates the origin is at y = 1 while the past light cone of the singularity is at y = 0. Though we 
have introduced the coordinates T and y for convenience, the criterion for smoothness is that A, B and <f> be smooth 
functions of u and w. The background solution is smooth provided that n is a positive integer. 

To study the global structure of the background spacetimes we make another coordinate change y = x r ' 
becomes 



The metric 



ds 2 = 
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(2ndx ~ x dT) dT + - (l - x 2n f dQ 2 
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The maximal extension of the spacetime is provided by this metric with — 00 < T < 00 and — 1 < x < 1. The Ricci 
scalar is proportional to e 2T , and so T — 00 is a curvature singularity. The regular center is at x = 1, and the past 
lightcone is at x = 0. Lines of constant x (trajectories of the homothetic vector field) are timelike for x > and 
spacelike for x < 0. The rings of constant x and T are closed trapped surfaces for x < 0, and so x = can be 
interpreted as an apparent horizon. 

Note that the x < region of the spacetime contains trapped surfaces, while a critical solution (since it lies on 
the boundary of those spacetimes with and those without trapped surfaces) cannot itself contain a trapped surface. 
Thus it is only the x > part of the spacetime that matches the corresponding part of the numerical critical solution. 
That part of the critical solution that is outside the past light cone of the singularity is not given by our background 
solution. 

In 2+1 dimensions, the cosmological constant is also necessary for black hole formation, and this goes beyond the 
approximation A = here (and beyond our current understanding). In a related complication, in the presence of a 
cosmological constant, a black hole will always eventually form and capture all the scalar field matter. The threshold 
behaviour and mass scaling studied by Pretorius and Choptuik applies only to prompt black hole formation by a 
scalar field pulse that has not yet reached the AdS timelike null- infinity. 
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III. PERTURBATIONS 



We now consider perturbations of the background solution. Here we use a 6 to denote a perturbed quantity. Any 
perturbation mode is of the form 

SB = e (k -^ T b{y), (10) 
6<t> = e kT H(y), (11) 
SA = e kT a(y), (12) 

where k is a constant. The perturbation grows as the singularity is approached if k > 0. 

These perturbations must satisfy the linearized versions of equations (j^-||) which we now solve in turn. Substituting 
the ansatz ( |l(i| ) into the perturbation of equation (||), we find 



yb" + (2k - 1)6' = 0. (13) 



The general solution is 



b = c Q + Cl (l-y 2 - 2k ), (14) 



for arbitrary constants cq and c\ . We may also ask which infinitesimal gauge transformations x^ — > x^ + £^ create 
metric perturbations of the form ( [l0| - |l2] ). It turns out that they form a 2-parameter family with the same parameters 
Co and C\. The perturbation b is therefore pure gauge. In order to study regularity at the origin, we introduce 
Cartesian coordinates t = (u — w 2n )/2 and f = (—it — w 2n )/2. The background metric then approaches the usual 
Minkowski form at the origin. All metric coefficients, and the scalar field <j), are regular at f = 0, and are even 
functions of f. The gauge transformation takes the form 

e = -i(co + Ci)H + f) l ~ k + Cl (-i - rf' k l (15) 
e = (c () + c 1 )(-t + f) 1 - k -c l (-t-f) 1 - k . (16) 

This is a regular vector field at r = if and only if Co = 0, and we assume this from now on. We are left with a family 
of linear perturbation gauges parameterized by c\. We shall later fix ci by further regularity considerations. 
The perturbation of equation (^|) now becomes 

\y(\ y 2 )H" + [k-(k + l)y 2 ]H' kyH = [k(l + y) - y + y^l - k(l + y))] . (17) 

1 (i + y) 

The solution of this equation that is regular at the origin (y = 1) is 

2cci 

H = -— i- (1 + y 1 - 2 ") + c 2 F(k, 1/2, 1, 1 - y 2 ), (18) 

l + y 

where C2 is a constant and F is a hypergeometric function. When C2 = 0, the perturbation simply results from 
applying an infinitesimal coordinate transformation to the background and is therefore pure gauge. 

We now address the issue of smoothness of the perturbation at y = 0, the past light cone of the singularity. We 
first consider the particular cases k — 1 and k ~ 1/2. In the case k = 1, we have F (1, 1/2, 1,1 — y 2 ) = 1/y. Therefore 
regularity at y = requires that H — for the case k = 1. In the case k = 1/2, we have F(l/2, 1/2, 1, 1 — y 2 ) — 
(2/ir)K(^l-y 2 ) where K is an elliptic integral. The quantity K(\J\ — y 2 ) diverges at y = 0. Since the first term on 
the right hand side of equation ( |l8| ) is regular at y = for k = 1/2, it follows that for k = 1/2 a regular perturbation 
must have C2 = 0, which is pure gauge. 

Now we consider general values of k. Using an identity for hypergeometric functions (section 9.13 of ||]) we have 

F(k, 1/2, 1, 1 - y 2 ) = J_ C2 { ~ k) k) F(k, l± + k,y 2 ) + ^Z^ y ^F(l y 2 ). (19) 

Here T(x) denotes the gamma function. This formula applies except in the case where k is half of an odd integer. In 
that case, one of the two r-functions has a pole. The correct formula then involves a In y term, and is therefore not 
regular at y = 0. All odd half- integer values of k must therefore be ruled out. 

The hypergeometric functions on the right hand side of equation ( [l9| ) are power series in y 2 and since y cx w n they 
are power series in w. However, due to the expression y 1 ~ 2k the right hand side of equation (091) will not consist of 
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integer powers of w unless k = m/(2n) for some integer m. We now consider the possible values of m. The cases 
to = n (k = 1/2), m = In (k = 1) and m = (2p — l)n for p integer (k odd half-integer) have already been ruled out. 
In considering the other cases, it is helpful to introduce the quantity C3 = C2T(fc — 5) / (^/tt T (k)) . The solution for 
H (x) near x = is then 
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For m < 77 these are all positive powers of X. For 77 < m < 2n., we have one negative power, but it is canceled when 
c\ = C3/ (2c). This means that there is only this one gauge choice Ci in which this perturbation is regular. For 777 > 277 
we have a negative power x 2n ~ m which cannot be canceled, and we must therefore have m < 2n. 

We now consider the perturbation of equation (0) which (after some straightforward but tedious algebra) becomes 



da 
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dy (1 + yf 
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Here F' is an abbreviation for (d/dy)F(k, 1/2, 1, 1 — y 2 ). 

Regularity of the perturbation imposes the condition that the quantity inside square brackets in equation j2l| ) 
vanishes at y = 0. Therefore, for each term we need only consider that part of it that does not automatically vanish 
in this limit. We first consider the case where m < n. Then, y 1_2fe vanishes at y = but y~ 2k does not. Expanding 
the quantity inside square brackets in equation (pl|) to the appropriate order we find 
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Here o{y) denotes a quantity that vanishes when y = 0. In general, the constant C3 can be chosen so that the quantity 
inside square brackets in equation ( p2] ) is 0(77). However, for 2nk = 1 (i.e. m = 1) the choice of c 3 that does this 
is C3 = and thus the perturbation is pure gauge. We then arrive at the result that a physical perturbation cannot 
have 777 = 1. 

Now consider the case 77 < m < 2?7. Then both y~ 2k and y 1_2k are nonvanishing as y —> 0. In addition, we have 
the restriction that C3 = 2cci . Then expanding the quantity inside square brackets in equation ( pl|) to the appropriate 
order we find 



da ci(k 

ny— = 

dy (1 + y) 



^{2nk + 1-2/7) [y~ 2k + 2i/- 2k + o(y)] 



(23) 



Thus, the only way to obtain regularity of the perturbation is to have 2?r,fc = 2n — 1 or in other words m = 2n — 1. 
Note that there is a mode with m = 2n — 1 only when this value is greater than 77, that is for 77 > 1 . 

Regularity at the origin y = 1 does not give rise to further regularity conditions beyond the ones already discussed. 



IV. CONCLUSIONS 



Perturbation modes depend on T as e kT . As a result of imposing regularity conditions, we have found the discrete 
spectrum of growing perturbations k = m/(2n) 1 where m is an integer with the following restrictions: 777 > 1 and 
either m = 2n — 1, or m < n. This means that the CSS background solution with index 77 has 77 — 1 unstable (k > 0) 
modes. In particular, the solution with n = 1 has no unstable modes. The 77 = 2 solution has one unstable mode 
with 777 = 3. The solutions with 77 > 2 all have one unstable mode with m = 2n — 1 and 77 — 2 additional unstable 
modes with 1 < m < n. 

We now compare our results with the numerical results of Ref. ||. The critical solution should be the one with 
one unstable mode. Thus our analysis indicates that the critical solution is the one with n = 2 and therefore that 
k = 3/4. In contrast, the analysis of curvature scaling in || indicates that 7 = 1.2 ± 0.05 which corresponds to 
k = 0.83 ± 0.04. Furthermore, direct comparison between numerical and analytic critical solutions indicates a much 
better fit with n = 4 than with 77 = 2. The 77 = 4 solution has an unstable mode with k = 7/8 = 0.875 which is in 
agreement with the k of Ref. However, our analysis indicates that it also has modes with k = 1/4 and k = 3/8 
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and therefore this solution cannot be the critical solution. Similarly, the n = 3 solution has an unstable mode with 
k = 5/6 which corresponds to a 7 of exactly 1.2 but it also has a unstable mode with k — 1/3 and therefore it seems 
that it cannot be the critical solution. 

Alternatively it may be that some additional condition needs to be imposed on the perturbations due to the fact 
that the past light cone of the singularity is an apparent horizon. Recall that the critical solution is only equal to our 
background solution within the past light cone of the singularity and must be something else outside the past light 
cone in order to avoid trapped surfaces. Thus this critical collapse situation differs from the usual one, and it may 
be that smoothness of the perturbation is not a sufficiently strong condition to impose. Such an additional condition 
might eliminate some of the modes that we have found. Some light could be shed on this issue by further numerical 
treatment of the original critical collapse problem. In particular, such a treatment could find the growing mode and 
compare it to the modes found in our perturbative treatment. 
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